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We study the planar-flow distributions of narrow, highly boosted, massive QCD 
jets. Using the factorization properties of QCD in the collinear limit, we compute 
the planar-flow jet function from the one-to-three splitting function at tree-level. We 
derive the leading-log behavior of the jet function analytically. We also compare our 
semi-analytic jet function with parton-shower predictions using various generators. 
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I. INTRODUCTION 

Observables sensitive to the substructure of energetic, ultra-massive jets hold great promise 
for distinguishing new physics signals from QCD backgrounds. Both the ATLAS and 
CMS experiments are pursuing studies relying on such observables in various new physics 
searches [i , (see also Refs. [3, 1], and Ref. ['] for an earlier search by the CDF collabora- 
tion). Boosted jets originating from electroweak gauge bosons [6], top quarks [7, 8], Higgs 
bosons [ ], and even new physics particles [10, 11] are all of interest as targets of searches 
by the Tevatron and the LHC experiments. It is therefore important to be able to dis- 
tinguish them from QCD jets. Recent reviews on substructure techniques, experimental 
status and new physics searches include Refs. [12-16] and references therein. 

One way to characterize jet substructure is to consider observables which are functions of 
the energy flow within the jet, namely the energy distribution as measured by the detector 
(see Ref. [17] for a recent systematic classiflcation) . In this work we consider the hadronic 
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collision 



(1) 



where H^, Hb are the initial hadrons, and J is a jet with momentum given by and rj, 
with size R determined by the jet algorithm, and characterized by an energy-flow observ- 
able O such as the jet mass m. Note that O can stand for multiple energy-flow observables. 
We focus on narrow, highly-boosted jets, and consider the inclusive differential QCD cross 
section for this process. 



Such cross sections are usually computed using parton-shower codes, which offer much 
less insight into results than analytic computations. They also require significant compu- 
tational resources. 

In this article we follow a different approach, based on the collinear factorization properties 
of amplitudes in perturbative, massless QCD. Factorization allows us to focus on a single 
jet, ignoring to leading order the rest of process, and to compute the differential jet 
substructure distribution semi-analytically. In some limiting cases, we can compute the 
distribution completely analytically. As we demonstrate below, this method becomes 
useful when analyzing jets with sufficiently high px {pt ^ 1 TeV), a window recently 
opened at the LHC (see for example Refs. [1, 2]). 

There are various ways to define jet shapes. In the context of new-physics searches, a 
particularly useful way to characterize jet shape and substructure observables is by the 
first non-trivial order at which they appear in fixed-order perturbation theory [17-19]. 
Consider, for instance, Higgs boson searches: at leading order the Higgs boson decays to 
two partons with large invariant mass. The relevant non-trivial substructure observables 
must distinguish two-prong jets from the broad spectrum of all jets. Given a Higgs-boson 
mass and pt, its decay kinematics are fully determined by one additional continuous 
variable, such as the ratio between the two decay-product momenta [ ] or the opening 



dp^ drj dO 
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angle between them [J.8]. The first non-trivial QCD background arises from corrections 
in which the jet is made up of two partons. If the Higgs boson recoils against a jet, this 
background first arises at next-to-leading order (NLO) in two-jet production, when 2—7-3 
processes are taken into account. This contribution is of 0{al) in fixed-order perturbation 
theory. 

In this example, the signal distribution at fixed is fully characterized by the jet mass 
and jet angularity [18, 20]. We may think of these two quantities as replacing partonic 
(and therefore unphysical) parameters with physical infrared- and coUinear- (IRC) safe 
jet-shape observables [18]. When working to leading order (LO) in the jet mass, the an- 
gularity distribution is the only independent jet-shape observable that can separate the 
signal from the background. The corresponding LO distribution, given a mass cut, can be 
computed analytically both for the signal and background [ 1 8], using the collinear approx- 
imation, which is adequate for narrow massive jets. The difference between the signal and 
background distributions turns out to be modest, because both the QCD and Higgs-boson 
angularity distributions are monotonically decreasing functions between identical limiting 
values of the angularity variable. The same result obtains when considering the ratio of 
momenta or other kinematical variables, because they are all fully correlated with the 
angularity distribution. The similarity of the bounds on both the signal and background 
does yield a sharp prediction of this apparently naive picture. This prediction was qual- 
itatively verified experimentally in the CDF collaboration's study [ ] of high-px massive 
jets. 

The next example of interest involves studies of high-px top quarks. At leading order, each 
top quark decays to three partons; in the decay, each parton pair typically has a large 
invariant mass. The same configuration is also relevant to studies of new physics [ , ], for 
example of gluino decay in i?-parity violating scenarios [ : ]. A useful jet-shape observable 
in such a study is the planar flow [18, 22, 23]. If we focus on studying one top quark 
out of the produced pair, the leading QCD background is two-jet production where we 
constrain one of the jets to have a significant planar flow. This background first arises at 
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next-to-next-to-leading order (NNLO) in two-jet production, when 2 — 4 processes are 
taken into account. The contribution is of 0{aj) in fixed-order perturbation theory. 

While the planar-fiow distribution of a top-quark jet at leading order can be computed 
straightforwardly from its known matrix element, the planar-fiow distribution of the QCD 
background has not yet been computed for narrow massive jets. The corresponding dis- 
tribution was presented in Ref. [ .] with a rather limited sample size due to the limited 
statistics of massive boosted jets at the Tevatron. LHC experiments have collected a 
much larger number of massive jets, which should allow a more precise measurement of 
the planar-fiow distribution. 

Our main purpose in this article is to compute the planar-fiow distribution of narrow 
massive jets to leading order in QCD. We use the coUinear approximation, in which we 
approximate the matrix element for such jets by 1 — t- 3 collinear splitting functions. 
Motivated by boosted-top studies, we take the jet mass to be roughly the top-quark mass. 
Our approximation is relevant when the jet px is substantially larger than the jet mass; 
we take the to be (9(1 TeV). Our computation assumes a jet algorithm that produces 
approximately circular jets with radius R in the pseudorapidity-azimuthal angle plane, 
but is otherwise general. The use of the collinear approximation also requires that the 
jet radius not be too large; we take R = 0.4. In parton-shower simulations to which we 
compare, we use the anti-fcx algorithm, with the same jet size. As already mentioned, this 
year's tt-resonance searches are already exploring this range of parameter space. 

Jet shape observables can be viewed as moments of the energy distribution within a 
jet [17]. They are highly susceptible to contamination from pile-up and other sources 
of soft radiation, especially for larger cone sizes [24, 25]. Such contamination is a major 
concern at present, with more than 20 interactions on average at each LHC bunch crossing. 
This number is expected to grow even larger in future runs. Various techniques [24, 
26] allow one to estimate and subtract pile-up contributions. Approaches in which jet- 
substructure analyses and searches are done in a way which is inherently less susceptible 
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to such contamination would offer desirable alternatives to contamination subtraction. 

Two main classes of alternative approaches have emerged: 'filtering' [9] (see also [27]) 
and 'template overlap' [ ]. In the former, a measured jet is declustered and its soft 
components are removed. This leaves only its hard components to be reclustered into the 
'filtered' reclustered jet. In the second approach measured jets are not manipulated, and 
are instead compared to a set of templates built according to a chosen (computed) fixed- 
order distribution of signal jets. The comparison makes use of an 'overlap function' which 
evaluates the degree of overlap between each measured jet and the set of templates. The 
reader will find a discussion of jet-substructure observables and experimental applications 
in Refs. [28-31]. 

For both alternative approaches, it would be useful to study distributions of the core 
(hard) component of jets. This is relatively straightforward for the signals, but much 
more challenging for the QCD background. The semi-analytic calculations we pursue here 
are a first step in this direction, as our results provide a semi-analytic understanding of 
the kinematical distributions of the hard component of massive jets with non-trivial three- 
body kinematics. For this purpose we also compare the result of our full 1 — 3 calculation 
with a similar calculation employing an iterated 1 — )■ 2 collinear splitting function in the 
approximation to the matrix element. We also discuss how various scale choices affect our 
result. Finally, we compare our results with parton-shower results, both with and without 
matching to tree-level matrix elements. We cross-check our calculations with a simple 
analytic expression for the planar-flow distribution in the small-planar-flow limit. 

The paper is organized as follows. In sections II and III, we define the planar-flow observ- 
able, introduce narrow jets and discuss various aspects of calculations of the jet function. 
In section IV, we relate the jet functions of narrow jets to collinear splitting functions 
for two-body jets. As an illustrative example we compute the leading-order jet function 
for the mass distribution using splitting functions. In section IV B we prove that spin 
correlations in the splitting function factorize in all cases that are of interest in this arti- 
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cle. In section V we use these methods to compute the leading-order combined mass and 
planar-flow jet function of three-body jets. For small values of the planar flow we obtain 
an analytic result below, while at large planar flow we rely on numerical integration. We 
also compare the jet function to one computed using an iterated 1 — )■ 2 approximation, 
discuss the behavior of the jet function at large planar flow, and discuss the sensitivity 
to scale choices. In section VI, we compare the semi-analytic jet function to parton- 
shower calculations with and without matching to fixed-order matrix elements. We also 
discuss hadronization corrections, and corrections from terms in matrix elements beyond 
the collinear approximation. We present our conclusions in section VII. Six appendices 
furnish a variety of technical details. 



We are interested in two energy flow observables — the jet mass and planar flow. The jet 
mass squared is. 



where i refers here to any parton (or hadron or tower or topocluster in a more realistic 
experimental context) inside the jet. For a hadronic collider we may define the planar flow 
as follows. Define first a 2 x 2 momentum shape tensor. 



where (Ar^j, A0j) = (r/j — r]^c\.^4'i ~ 0jct) are the pseudorapidity and azimuthal-angle dif- 
ference of each jet constituent from the jet axis. We take all constituents to be massless. 
This form is manifestly boost-invariant for boosts along the beam axis. 

The planar flow is then defined by 



II. PLANAR FLOW 




(3) 





Pf = 



4detX 



(5) 
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One can easily verify that < Pf < 1, and that it vanishes for two-body jets, receiving 
its leading contribution from three-body jets.^ The latter property means that the planar 
flow is potentially useful for distinguishing the QCD background from top jets. The value 
Pf = arises when the partons lie on a line in the detector plane. In particular, it 
will vanish for three-body jets when a parton becomes soft, or when two partons become 
collinear. Three-parton configurations that are symmetric about the jet axis have Pf = 1. 
Configurations with n partons symmetric under rotations by 27r/n radians around the jet 
axis also have Pf = 1. 

For the sake of convenience, in the theoretical calculations that follow we will focus on 
central jets (that is, with very small pseudorapidity) . For such jets we can work in terms 
of the angles 9 and (being the polar coordinates around the jet axis), set ?7 = 0, and 
exchange pt for P (the overall jet momentum). For narrow, central jets, we can rewrite 
the momentum-inertia tensor as follows, 

M ^ y^ P±^P±^ , k,l = 1,2, (6) 

iSjet 

where p±^i^k is the kth component of the transverse momentum of constituent i with respect 
to the jet axis. 

As we review in Sect. IV, the leading-log behavior of the jet function for the jet mass is 
given by 

Jf{m';p^,R) ^ ^^log , (7) 

where in this case Cg = Ca = 3 and Cq = Cp = 4/3, and the jet radius R = a/ A?]^ + A^^. 
In this article we compute the mass and planar- flow jet function at leading order in a^, 
for narrow QCD jets. In the limit of small Pf we will obtain an analytic result for the 



^ Notice that because TrI is proportional to the jet mass [ ' ], the planar flow is only well defined for 
massive jets. 
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leading-log behavior of the jet function. For small m/px-R it is given by the expression 

c, = 2Cl, c, = Cf{Ca + Cf). (9) 

Away from the limit of small planar flow, we will compute the planar-fiow jet function 
semi-analytically, using numerical integration to obtain the final result. We will show 
that there is a physically interesting regime of parameters, with the jet mass near the 
top-quark mass and with 0.4 ^ Pf ^ 0.95, in which our result has rough agreement with 
the parton-shower simulations. (As shown below, the parton-shower results of different 
tools do not fully agree with each other, which however is not the focus of this study.) 
We expect our results to be useful for understanding how to refine methods to distinguish 
highly-boosted top jets from the QCD background at the LHC. 

The planar-flow distribution was measured by CDF [ :] and ATLAS [ >], but with large 
statistical uncertainties, and using too large a cone size and too big a mass-to-momentum 
ratio to be compared with our results. At the LHC, accurate measurements of planar- 
fiow distributions are difficult due to pile-up effects; but we may expect them to improve 
significantly over time. In principle one can 'refine' jets in a controlled manner (by applying 
filtering; by using the template overlap method and then looking at the parton distribution 
of the peak templates; by looking at events with a small number of vertices; or by using 
other methods for pileup subtraction), and thereby isolate the hard part of the measured 
jet in order to compare with theoretical predictions even in the presence of incoherent soft 
radiation. 

III. NARROW MASSIVE JETS 

Let us consider jets at the parton level. If we take a jet to be narrow {E? <^ 1), the 
partons will be approximately coUinear. As we review in the next section, when massless 
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partons become coUinear the QCD cross section factorizes partially, and we can write 
(schematically for now), 



In the cross section on the right-hand side, the partons making up the jet are replaced by a 
parent parton of type /, which can be either a gluon or one of the massless quarks. The jet 
substructure is encoded in the jet function Jj, which has a simple physical interpretation: 
it is the probability distribution for the parent parton of type / to evolve into a jet of 



Throughout our study we will be agnostic about the specific jet algorithm used in the 
analysis, and will assume only that it produces approximately circular jets with radius R 
in the ■r]-(j) plane. 

As we review in the next section, at fixed order in as the jet function can be computed 
from splitting functions [33], universal functions that govern the behavior of the squared 
matrix element in the collinear limit. In this limit, the squared matrix element factors 
into a product of a splitting function and a squared matrix element with lower multiplic- 
ity. In general, the factorization is not complete, due to the dependence of the splitting 
functions on the spin of the parent parton. We will show, however, that for all energy-flow 
observables the spin dependence does factorize. 

The fixed-order splitting function is singular in the limit where partons become soft or 
collinear. In Eq. (10) this singularity appears as a divergence of the jet function in the 
infrared limit of the observable O (for example taking m — )■ in a two-body jet). Resum- 
ming higher-order (perturbative) corrections cures the divergence, and replaces it with 
a peak at a finite value of O (for reviews see for example Refs. [34-37]). A fixed-order 
calculation is accordingly unreliable when we get close to the infrared limit. This problem 
can be avoided by considering only values of the observable O that are far away from 
the infrared limit, compared with the peak position in its distribution. In particular, we 
will always take the jet mass to be 'large enough' in this sense. The peak in the jet-mass 




(10) 



size R that has O = Oq. Accordingly, its integral is normalized to unity, / dOJf{0) = 1. 
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distribution, mpeak, is roughly near its average at m^ve = cxsPiR [12], and we will take 
our jet mass to be much larger, m ^ ^pcakj where we expect the fixed-order perturbative 
calculation to be reliable. 

Beyond higher-order perturbative corrections requiring resummation, there are also non- 
perturbative corrections (in the form of hadronization), which become important in the 
infrared and tend to smear the jet function. If some of the fixed-order partons have 
transverse momentum relative to the jet that is small compared to the characteristic 
transverse momentum of the smearing effect, the final jet function will be dominated 
by the latter effects. Keeping the observable O away from its infrared limit avoids this 
problem as well. On the other hand, in order for the coUinear approximation to hold, 
we cannot stray too far from the infrared limit. Finally, in order for the approximated 
distribution (in the coUinear limit and at fixed order in a^) to be valid, we should not get 
too close to kinematic boundaries. 

For our approximation to be reliable, we need a range of values for O that obeys these 
constraints. Applying the constraints to the jet-mass observable, and requiring that we 
have a non-empty range of validity for the approximation, necessitates considering jets 
with sufficiently high px- We therefore consider only highly-boosted jets. For a general 
observable, the existence of such a range of validity is less clear. We will show later that 
to reasonable accuracy, a non-trivial range of validity does indeed exist for the mass and 
planar-fiow jet functions with collision parameters typical of the LHC. 

IV. QCD JET-MASS DISTRIBUTION IN THE COLLINEAR LIMIT 

In this section we compute the jet-mass distribution for massless QCD in the coUinear 
limit using the splitting function [33], for both quark and gluon jets. Consider again the 
hadronic collision. 



HA{qA) + HB{qB) ^ JiO-p^.T]- R) + X , 



(11) 



13 



where Ha,Hb are the initial hadrons with momenta qA,qB, and J is a jet of cone size 
R and given and 77. The jet is further characterized by an energy-flow observable O. 
For simplicity, in this section we will take O to be the jet mass squared, m^; in the next 
section we will consider also the planar flow. 

The factorized cross section is given by, 

where (pf^,4>fB parton distribution functions. For narrow jets daj^fg^jx further 

factorizes at leading order into a jet function, times a cross section in which the jet is 
replaced by a single parent parton. 



dpi^d^'dm'^ (a^^. xb, Pt, V, R) ~ ^ '^^dp^d^^ ^b,Pt, v) Jfirn;pT,r];R) . (13) 



The sum is over the type / of the parent parton, which can be a gluon or a quark with 
a specific flavor. The relation (13) is due to the factorization of the QCD matrix element 
in the limit where two or more partons become collinear. The universal function in the 
factorization is proportional to the splitting function. We therefore seek to express the jet 
function in terms of the splitting function. 



A. Quark Jets 

Let us first consider the case where the parent parton is a quark. This case is simpler 
because the quark splitting function does not depend on the helicity of the parent parton. 
Here the cross section factorizes completely in the collinear limit. Consider the matrix 
element for the scattering of two into n massless QCD partons. 

Here pi denote parton momenta; the fi, parton types; the q, their colors; and the Sj, 
their helicities. Outgoing particles are indexed by z = 1, . . . , ra, and incoming particles by 
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i = n + l,n + 2. Define the abbreviation |A^2^>nP = l-^/i,...,/„+2(Pi5 • • • )Pn+2)P for the 
squared matrix element summed over cofor and hehcities (averaged in the case of incoming 
particles — we will leave this distinction implicit from now on). 

Consider the limit in which two outgoing partons, say pi and p2, become collinear. The 
leading contribution to |A^2^nP in this limit is from diagrams in which the two outgoing 
partons originate from a single parent parton with momentum p = Pi + P2, and with a 
type / that is uniquely determined by the splitting process / — /i /2- In the collinear 
limit the parent goes on shell, leading to a l/p"^ singularity. When the parent is a quark, 
the squared matrix element f actor izes as we approach the limit, 

S12 

(15) 

where S12 = (pi + ^2)^, and -P/1/2 is the spin-averaged splitting function [33, 38], given 
in appendix A. In the squared matrix element on the right-hand side, the two collinear 
partons are replaced by their parent parton. For a gluon jet, the splitting function depends 
on the helicity of the parent parton, and the factorization is not as simple as in Eq. (15); 
we consider this case in the next subsection. 

The fixed-order differential cross section is given in terms of the squared matrix ele- 
ment [39], 

1 d^Pi 



1 



n+l-C/n+2 .^-^ 



2Ei (27r) 



\M2^nr S%Pi,tot - Pi,tot) . (16) 



Using Eq. (15) and making a change of variables (pi,P2) — ^ {Pi^P = Pi + P2) it is easy to 
see that near the collinear limit, 



47ra, E d^pi 
da2^n ^ da2^n-i X - — ————^Pf,f,{pi,p2] 

S12 rjirj2 [iTiy^ 



(17) 

P2=P-Pl 



where d(j2^n-i is defined in terms of the matrix element on the right-hand side of Eq. (15). 
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Comparing with the postulated relation (13), we may now write down the jet function, 

X 6{m\pi,p2) - m') e(i? - ^i) Q{R- 62) ■ (18) 

Here, 9i are the angles of the momenta pi with respect to the jet axis p, and the step 
functions Q{R — 6i) are put in by hand to enforce^ a cone of size R. The sum is over 
allowed splitting processes, and the symmetry factor Sf-^f^ corrects the over-counting of 
identical parton configurations in Eq. (18): it is 2 when fi = /2, and 1 when they are 
different. We stress that Eq. (18) is valid only to leading order in as at the partonic level. 
The full jet function receives corrections at higher orders in perturbation theory (some of 
which require resummation) as well as from non-perturbative effects. 

It is now a straightforward exercise to substitute the quark splitting function into Eq. (18) 
and compute the jet function, assuming _R ^ 1. The quark splitting function is [33, 38] 

Pq^,q{z) = Cf ^ ^ , (19) 

where Cp = {N^ — l)/2Nc, and z = Eg/E is the emitted gluon's energy fraction, in our 
approximation. Using Eq. (19) and solving the integral in Eq. (18), we find the leading-log 
expression for the quark jet function, valid for m ^ pR 

Um^-p-R)^^\\oJ^^) . (20) 

This result was derived in Ref. [22] using slightly different terminology but a similar limit. 
As we explain in detail in appendix F, the form of Eq. (20) can be alternatively obtained 
by rewriting the mass as a function of the emission angle and z, and replacing z with 
m as the integration variable. This makes explicit the overall l/m? dependence, and the 
integration over the angle within the allowed kinematical boundaries further leads to the 
log in Eq. (20). 



^ This procedure is expected to be compatible, up to higher-order corrections in R, with any jet algorithm 
that produces approximately circular jets [18]. 
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When m — )■ 0, the jet function (20) diverges. This is an infrared divergence, resulting from 
partons becoming soft and colhnear. In this hmit higher order contributions in «<, become 
important, and after resumming them the singularity is exponentially suppressed. Thus, 
the full jet function vanishes in the massless jet limit (see for example Refs. [34-37]); and 
because it decays at large mass, a peak of the jet mass is expected to arise at low jet mass. 
The result in Eq. (20) is therefore only reliable when mpeak '^ttl <^ pR. The divergence in 
Eq. (20) also renders the distribution non-normalizable. Nevertheless, in the regime where 
our approximation is valid we expect our result to match the full jet function including 
its overall normalization. In other words, in this regime it gives the (leading) probability 
distribution of a quark to evolve into a radius- i? jet with mass m [18, 22]. This behavior 
was verified experimentally, at least qualitatively, in Ref. [5]. 



B. Gluon Jets 



We now turn to the case of gluon jets. Unlike the quark case, the gluon splitting functions 
have a non-trivial dependence on the helicity of the parent parton. As a result, the colhnear 
factorization of the squared matrix element is incomplete. For the computation of the jet- 
mass distribution the incomplete factorization does not alter the leading-log result: in this 
case, the spin correlations are nonsingular. However, we would like to go beyond this and 
establish a general result that will be useful in the rest of our study. The spin correlations 
always factorize when one considers distributions of energy-flow observables deflned for a 
single jet. We will prove this for the 1 — 2 splitting function, and the proof can be easily 
generalized to the 1 — )■ 3 splitting function. 

In the presence of spin correlations, the colhnear factorization (15) of the squared matrix 
element is no longer correct, and instead we have [iS], 

ss' 

(21) 
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where now P is the hehcity-dependent sphtting function, and T is defined in terms of the 
matrix element M.2^n-i by 



-^/i,...,/n+i [Pi,- ■ ■ ,Pn+lj J^f^,...J„+, [Pl,---,Pn+l 



S2,...,Sn+l 
Cl,C2,...,Cn+l 



(22) 

T is essentially the matrix element A^2^n-i, squared, except that there is no sum over 
the helicity of the jet's parent parton. 

The differential cross section no longer factorizes as it did in Eq. (17). However, in the 
jet-function definition the phase-space integral includes an azimuthal integral around the 
jet axis. For example, changing the integration variables in Eq. (18) to polar coordinates 
{pi,6i,(f)i), i = 1,2 (relative to the jet axis p), the integration over rotations around the 
jet axis is described by the variable (p = (pi, keeping (<^2 — 4>i) fixed. Therefore, for any 
observable that is invariant under rotations around the jet axis the integral picks out the 
part of the splitting function that is invariant under such rotations, which is precisely the 
spin-averaged splitting function. This includes all energy-fiow observables defined in terms 
of a single jet, because for such observables there is no preferred direction in the detector 
plane that can break the rotational symmetry. We may therefore replace P^^' — ?■ 5^^' P, 
where P is the spin-averaged splitting function. Noting that 'j''^'^' = |M2^n-iP, the 
rest of the computation follows through as in the previous section, and we conclude that 
the jet function for gluon jets is given by Eq. (18), just as for quark jets. 

We conclude this section by computing the leading-log part of the gluon jet-mass function. 
For <^ 1 and ^ {pRY, the leading contribution comes only from the g ^ gg 
splitting function. 



Pgg{z) = 2Ca 



z 1 - z , 

+ + z{l-z) 

1 — z z 



(23) 
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where Ca = Nc, and we find 

J,(m^;p;i?)^^i,logf^) , (24) 
again in agreement with Ref. [22]. 

V. PLANAR-FLOW JET FUNCTION 

In this section we compute the planar-fiow jet function J/(m^, Pf; px; -R), which receives 
its leading contribution from three-body jets. This jet function factorizes from the rest of 
the cross section when we take the "triple" coUinear limit, in which three partons become 
collinear simultaneously. The limit is analogous to the ("double") collinear limit that we 
studied in the last section. 

We consider 2 — )■ n scattering with matrix element M., where three outgoing massless 
partons with four- momenta Pi,P2,P3 are collinear. The parton energies are denoted by 
Ei. The jet made out of these three partons has energy E = Yl^=i momentum p and 
mass m. Finally, we define the energy fractions Zi = Ei/E, and also, 

Sij = (pi + Pjf = 2pi ■ pj , z, j = 1, 2, 3 , 
Si23 = (Pi + P2 + PsY = m"^ ■ (25) 

In the collinear limit, the squared matrix element of the scattering process factorizes as 
follows [38], 

l-^/i,/2,/3,...,/n+2(Pl,P2,P3, • • • ,Pn+2)P ^ "' ^!l...,/n+2(P'P4, • • • , Pn+2) Pfff, f,ipi , P2 , Ps) ■ 

^123 

(26) 

Here, Pfff^f^ is the one-to-three sphtting function [38, 40] for the splitting / — )■ /1/2/3; the 
type / of the parent parton is determined by fiavor conservation. The splitting function 
depends on the outgoing momenta and each parton /j is associated to a momentum p,. 



19 



As in Sect. IV B, performing the azimuthal-angle integrals will make spin correlations 
disappear in the jet function. We may thus replace 'T'^^' Pjff^j^ by |M2^„_ipP/i/2/3, where 
Pfif2f-i is the spin-averaged splitting function. In the jet function we will need to sum 
over all possible final partons and, as before, take into account the symmetry factor that 
compensates for an over-counting of identical-parton configurations. Let us define then 

Pf= X] Q Phf2h ' (27) 

where Sf-^j^f^ equals n\ for a parton configuration that includes n identical partons. Ac- 
cordingly, 

Pg = ^^999 + ^fPgqq ) (28) 

Pq = ^P99q + l^ Pqm + {Nf - 1)Pqq, , (29) 

where Q denotes a quark with a different flavor than that of the parent q. The spin- 
averaged splitting functions Pf^f^f^ are given in appendix B. 

The jet function is given by a straightforward generalization of Eq. (18), 

Pf{pi)5[m^{pi) - m^) 

X <5(Pf(p-;) - Pf) 5'{p- ^p,) . (30) 

As discussed in detail further below (see subsection VE), the two strong-coupling factors 
are evaluated at different scales^: we take the first scale to be m, corresponding to the 
first 1 — )■ 2 splitting, and we choose the second scale to approximate that of the second 
splitting according to an ordering scheme defined below (obviously for "Mercedes" -like 
configurations the ordering does not make a difference). 



J/(m^Pf;p;i?) 



n 



-1 



^ For the sake of brevity we allow ourselves to be sloppy in our notation: generally ^ in Eq. (30) depends 
on the parton configuration in the integrand, so as(/Li) should really appear inside the integral. 
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An equivalent definition of tlie jet function (30) is provided by averaging tlie splitting 
function over parton permutations, 

wliile restricting tlie integration domain to 

/oo /• 
d^Pid^P2d^P3 — > 3 / d^pid^p2d^p3 . (32) 

■oo JEi<E2,E3 

The motivation for this alternative formulation will become clear below, but the basic 
idea will be to keep one parton (say, ps) away from its soft limit — )■ 0. We can do this 
without loss of generality as our jet is assumed to be of a large mass. 

Let us now integrate over ps in Eq. (30) by using momentum conservation. In the re- 
maining integral let us switch to spherical coordinates relative to the jet axis p, {pi, 9i, (pi), 
z = 1,2. 

Next, we extract the leading-order expressions for the integrand in the narrow-jet expan- 
sion 9i < R 1. The resulting expressions appear in appendix C. Using these, as well 
as the explicit splitting functions, one may easily check that the integrand in Eq. (30) 
depends on the azimuthal angles (pi only through the combination cos((/)i — (p2)- We make 
the change of variables J^"" d(pid(p2f{cos{(pi — 02)) = 47r d(pf{cos{(p)). 

We can now write the jet function as follows, 

Jf ^ ^ rde,de2 rd(P [ dp, dp2 ^-^^^Pm 5{m'ip^ - m') SiFm - Pf) , 

Jo Jo Jc PS 

(33) 

where the domain of integration is, 

C = {{pi,p2) \ 0<pi<p2,P3, 93<R} . (34) 



Here, ps and 63 are given (in terms of our integration variables) by eqs. (C5) and (C6) 
respectively. It can be easily seen that the resulting collinear expansion of ps, Eq. (C5), 
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breaks down for small ps. In other words, the soft limit of ps and the coUinear limit do 
not commute. Therefore, one cannot use the expansion (C5) in the vicinity of small p^. 
We avoid this complication by formulating the jet function in a way that avoids that part 
of the integration region, as in Eq. (32). 

Using Eq. (C6), the inequality 9^ < R is equivalent at leading order to 



Notice that we need m < pR to have a non-vanishing domain, in accordance with the 
narrow jet approximation. Very close to the kinematic boundary the leading contribution 
will be dominated by higher order terms; consequently, we must require m to be safely 
away from the kinematic boundary, ^ {pRY- 

Before proceeding, let us discuss the expected range of validity of our theoretical compu- 
tation. Take pt — P = 1 TeV. The collinear approximation requires that i?^ ^ 1; we 
choose i? = 0.4 for our analyses, the same as the smallest cone size used at the LHC for 
new-physics searches. As for the jet mass, we would like to stay well below the kinematic 
boundary at pR, but also well above the peak in the mass distribution (very roughly 
at aspR), so that higher-order effects in remain small. We choose m ^ 180 GeV, 
which satisfies both of these constraints. This choice is also well-motivated physically, as 
it is close to the top-quark mass. The singularity at small Pf would be eliminated by 
resummation of higher powers of log Pf; we would expect such a resummation to lead to 
a Sudakov-like exponential damping term, ~ exp (— log^Pf), similar to the case of the 
thrust or jet-mass distribution [ ]. To the best of our knowledge no computation of this 
effect has been carried out to date. In any case the focus of our present investigation is on 
massive jets with sizeable planar flow where higher-order corrections, resummed or not, 
have a subdominant impact. Qualitatively we expect the Pf distribution to be similar to 
that of the jet mass, namely, vanishing for Pf = 0, peaking at a small value of Pf, and 
falling gradually beyond that point. As we shall see in Sect. VI, the jet parameters chosen 
above lead to a Sudakov peak in the planar- flow distribution near Pf = 0.1 . We will 




) 



p,(R^ - el) - P2{R' - 91 



)>0. 



(35) 
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not restrict the range of planar-flow values in our discussion, but remind the reader that 
physically reliable results are expected within our framework only for planar-fiow values 
well above this peak. 

A. Analytic Leading-Log Behavior 

In this section we compute the jet function analytically in the limit of small planar flow 
and at fixed jet mass. For simplicity, we will take the second running scale /x in Eq. (30) 
equal to the first, m. We will find that the leading-log result, in the limit of small Pf and 
small m/pR, is given by Eq. (8). 

At Pf = 0, the integral in Eq. (33) diverges because the splitting function is singular. The 
singularity arises in regions of integration where a parton becomes soft, or two partons 
become coUinear (with respect to the third parton). In fact, the leading singularity in 
the integrand arises in the combined soft-collinear limit, where a single parton becomes 
both soft and collinear with another parton. In the splitting functions (see appendix B), 
the terms responsible for this leading singularity are those proportional^ to {sijZi)~^ or 

(SjjSjfc) 

At small (but nonzero) Pf, the leading contribution to the (finite) integral will thus come 
from the soft-collinear regions, which are disconnected in the domain of integration. In 
Eq. (33) there are two such regions, where parton 1 is soft, and collinear with either parton 
2 or 3. Due to the symmetry of the partons, it suffices to compute the contribution from 
one region. We will compute the contribution from the region in which parton 1 is soft 
{pi — 7- 0), and collinear with parton 2 {612 — )■ 0). 

One might have expected terms proportional to to lead to an even higher singularity. A careful 
examination of these terms, though, reveals that they are in fact less singular than the other ones 
mentioned above. 
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Let us choose more convenient variables in Eq. (33) to describe the colhnear hmit, 

9i = 9 , 02 = 0(1 + r cosco) , (f) = rsmu. (36) 



To leading order in r, we see from Eq. (C7) that 612 = Or, and the colhnear limit is now 
parameterized by r — )■ 0. 

In the soft-collinear limit {pi,r — )■ 0) the splitting function can be written as 

5 Ff{9;p,m) 



+ 



(37) 



where the dots include less singular terms; this is the leading singularity to which we 
alluded above. The functions Fj are given by 

CFTRNfW?p^e^ (m^ + p^e^) + Cl (m^ + mV^^ + P^O^f 



3pW^m^+pWy 
2CFm^ + {Ca + 2CF)mYe^ + {2Ca + CpWp^e^ + 2Ca/^^ 



(3^ 



(39) 



As we seek the leading Pf — )■ singularity, we will include only these terms in the splitting 
function. 

Using Eqs. (C5), (CIO) and (Cll), at leading order in pi,r and in the colhnear approxi- 
mation we have 



m 



^p-^p^p^O^r"^ sir? (jj 



Pf 



p- P2 ■m^{p - P2) 

The integration domain (34) becomes^ 



P3=P-P2 



C ^ < (Pl,P2) 



Pi > , (i?' - 9^)p2 <piR 



m 



(40) 



(41) 



The jet function (33) can now be written as 

Qa^p 



J 



f 



3™4 



d9 I dco I dpidp2 I dr x 



6 



P-P2 







m 



6 



c JO P-P2 

Ap^PiP26*r'^ siv? oj 



m'^{p - P2) 



-Pf 



(42) 



As we are interested only in the region where pi is soft, pi < P2,P3 is automaticahy satisfied at a fixed 
mass. 
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with a factor of 2 from the sum over the soft-coUinear regions. The upper integration hmit 
for r is a function of the other variables whose precise form is irrelevant to the leading 
behavior we are trying to compute. 

In the limit of small planar flow, the dominant contribution to the jet function is from 

the integration region near the soft-collinear singularity. We may therefore restrict our 
domain of integration to a small region around the singularity, defined by 

< pi < , < r < r"^^ , (43) 

where now both r™^ and are chosen constant and small. We expect that, in the 
small Pf limit, the jet function will not depend on these parameters. Indeed, we will see 
that r™^ and will drop out of the result. 

We can now solve the delta functions for pi and p2, 

m^Pf rn^p 
ApU^r^ sm uj + p^t^^ 

Taking into account the Jacobian and the appropriate domain, we find, 

n^m^PfJmVp^R {m^+P^O^yJo J^^^^r 

where the lower integration limit on r emerges from the previous upper integration limit 
on pi. Proceeding with the r integration, we finally find the leading- log approximation for 
the jet function. 



where 

..2^2 rR 9^Ff{0,p,m) 

As anticipated, the arbitrary integration limits r™^ and p™'^'^ do not appear in the leading 
term. When we include the first correction, the jet function takes the form. 



•^/^pfloglpf )+•••■ (48) 
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Here we will not compute Bf, which does depend on the values of r"^^^ and p'^^^. The 
remaining integral in Eq. (47) can be performed analytically, and the full result is given 
in appendix D. Here we record only the result to leading order in m/pR, 



where Cg = 2C^, Cg = Cp{Ca + Cp). This is the same leading-log behavior that we found 
in the case of the jet-mass function. 

For the relevant range of parameters as taken earlier, the contributions subleading in 
m/pR cannot be neglected. Taking p = 1 TeV, m = 180 GeV, R = 0.4, and evaluating 
the coupling at the jet mass scale, we find Ag = 0.683 TeV~^, Ag = 0.227 TeV^^. In 
the next section we will compare this result with a full semi-analytic evaluation of the jet 
function in the collinear approximation. 

In appendix F we show how, in parallel to the evaluation of the jet-mass distribution, the 
most singular term in the jet-Pf distribution can be obtained simply by iterating the 1 — )■ 2 
splitting function. The 1 /Pf factor arises from changing variables for the second emission- 
energy fraction while the log(Pf) factor results from integrating the angular variable within 
the kinematic constraints. 



Let us now return to the expression (33) for the jet function, and proceed without making 
further assumptions. We will carry out two integrations using the delta function con- 
straints, and compute the remaining integrals numerically to obtain the jet function. This 
computation is valid in the collinear approximation and for any planar flow. 

Let us perform the integration in Eq. (33) over pi and p2 by using the mass and planar-flow 
constraints. Solving these constraints, one finds two solution branches pf that should be 




(49) 



B. Semi- Analytic Evaluation 
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summed over, 

(±),. . _ 2mVg^gisin^(/) + m^pPfg2(g2-gicos0)±v/A 

p(^)(^l,^2,0)=pS^^(^2,^l,0), (51) 

where 

A = mYdldl sm20 {2p2^i02 [p^^i^2(l - cos(20) - 2Pf) - 2m2pf cos0] - m^Pf^} . (52) 
The jet function becomes 

where = d{p['\pi'^)/d{Ff, m?) is the Jacobian, and the new integration domain 

(including a reahty condition for pf^"^) is 

c:={(0i,02,</') |a>o, o<pp<p(^\p3(pf\^^,0), 

0<ei,2<i?, O<0<7r, e3(pf\^^,0) <i?}. (54) 

We compute the integral in Eq. (53) numerically.^ For the comparison with the analytic 
result we take both couplings at the jet-mass scale; when comparing with the parton- 
shower simulations we take the first coupling at the jet-mass and the second coupling at 
the dipole scale (see subsection VE and appendix E). Fig. 1 shows the results, including 
a fit to the predicted leading-log behavior. Table I compares the leading-log coefficients 
predicted in Sect. VA with those extracted from the fit. At planar-flow values below 
10~2 we find good agreement with the predicted leading-log behavior. At larger values 
of the planar flow the subleading behavior becomes important, and the fit must include 
the subleading coefficient Bf as well. We stress that in the small-Pf region, higher-order 



^ We obtained all semi-analytic results in this paper using Mathematica with adaptive Monte-Carlo 
integration. 
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effects are very important, so that the planar-ffow jet function obtained here requires 
resummation and cannot be compared usefully with experimental data. Above Pf = 10~^ 
the leading coefficient Af no longer agrees with the analytic prediction, which assumes 
Pf ^ 1, but the jet function does match the general form in Eq. (48). In this region, we 
must integrate the splitting functions numerically to obtain a semi-analytic prediction. 




(a) Small Pf (b) Large Pf 

FIG. 1: Comparison of semi-analytic jet functions (dots) with the purely analytic leading-log 
form of Eq. (48) (solid lines) and with a fit to the leading- log form (dashed lines). The quark 
functions are the lower (black) curves and the gluon ones are the upper (green) curves. Two 
ranges of the planar-flow observable are shown: (a) 10~'^ < Pf < 10~^ and (b) 0.1 < Pf < 0.8. 
The curve parameters are listed in table I; the purely analytic parameters are given in column 1, 
the fit parameters for (a) in column 2, and the fit parameters for (b) in column 3. In the purely 
analytic form, Af is given by Eq. (49) , and Bf = 1. In the fit to the leading-log form, Af is 
determined by fitting in both (a) and (b), while Bf = 1 is set to 1 in (a) and is determined by 
fitting in (b). The jet momentum is p = 1 TeV, with mass m = 180 GeV and size R = 0.4. 
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r\iifiiy lie 




Til Oil IJcLl^e -Ti 




0.683 


0.651 


0.413 




0.227 


0.205 


0.138 




1 


fixed to 1 


5.06 




1 


fixed to 1 


2.85 



TABLE I: Fit coefficients for the leading- log form in Eq. (48). The coefficients Aj are in units of 
TeV~^, where the Bf are dimensionless. The analytic result is only for the leading-order term 
(in small Pf), so that effectively Bf = 1. 

C. Comparison with Iterated 1 — 2 Splittings 

In this section we obtain a different semi-analytic approximation to the jet function, by 
approximating each 1 — 3 spHtting function by an iteration of two 1 — 2 splitting 
functions. For simplicity, we ignore spin correlations, which do not contribute at leading 
order in small Pf.^ This corresponds to the strongly ordered limit for the two splittings. 

Consider again the squared matrix element 

In the limit where pi and p2 are coUinear, that is have relative transverse momentum small 
compared to all other parton pairs, the matrix element (55) factorizes [cf. Eq. (15)] as, 

|M2-,„|^ - \Mf,j^^,„j„^^{p',P3, . . . ,Pn+2)\^ ■ ^^^/i/2(Pl,P2) , (56) 

where p' denotes the parent parton of partons 1 and 2, namely p' = pi + P2, and /' is 
determined by /i and /2. The splitting functions -P/1/2 are given in appendix A. Next, 
assume that p' and ^3 are also coUinear, with relative transverse momentum small com- 
pared to that of remaining parton pairs (but large compared to that of pi and P2). We 



^ The 1^-3 and (1 2)^ splitting functions can be shown to have the same soft-coUinear leading 
singularity. 
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then have the further factorization, 



\M, 



2-5>n 



\Mf 

(p,P4,...,Pn+2)h 



647r2a2 



■Pfif2{pi,P2)Prh{p,P3) 



(57) 



■S12S123 

where p and / are the jet's momentum and type, respectively. Comparing with Eq. (26) 
we can read off the iterated 1 — )> 2 approximation to the function Pf defined in Eq. (27), 

1 

- > - 

S12 



^hf2 



Pf'hi.P\P3)Phf2i.PuP2) ■ 



(58) 



Exphcitly, 



p(1^2)2 ^ ^123 

9 „ 

Sl2 



p(1^2)2 ^ ■Sl23 
Sl2 



^Pgg{p',P3)Pgg{Pl,P2) + A^/^gg (p', ^3)^55(^1, P2) 



+ 2NjP,gip\ps)Pg,{p,,P2) 



Pqg{p\P3)Pqg{Pl,P2) + ^Pgq{p' , P3)Pgg{Pl, P2) 



+ NfPgg{p',P3)Pggip^,P2] 



(59) 



(60) 



The strongly ordered approximation to the planar-fiow jet function is given by Eq. (53), 
with Pf replaced by Pj (symmetrized over parton permutations — see Eq. (31)). As 
in the case of the 1 3 splitting-function approximation to the planar-flow jet function, 
we may allow the argument scale of as to depend on the partonic configuration. We 
show the ratio of the jet functions in the strongly ordered approximation to those in the 
basic approximation of section VB in Fig. 2, where the scales for all factors of the strong 
coupling are fixed to the jet mass m. 
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(a) (b) 

FIG. 2: The ratio of semi-analytic jet functions evaluated with a splitting kernel approximated 
by an iterated 1—7-2 splittings, to those with the original 1 — )• 3 kernel. The parameters are 
p = 1 TeV, m = 180 GeV and R = 0.4, and the couplings are evaluated at the jet-mass scale. 



The strongly ordered approximation giving rise to iterated 1 — )■ 2 splitting functions 
implies a large hierarchy between the two splittings. This means it is valid in the limit of 
small Pf (keeping the mass fixed). We would thus expect the two jet functions to coincide 
in the small-Pf limit. In Fig. 2, the ratio does get closer to unity at small planar flow, 
but a gap appears to remain. This gap results from the discrete number of points chosen 
and the use of a linear scale; we have verified that it closes up when going further to very 
small Pf. In any event, a fixed-order calculation is not valid in this region; we should 
focus on the region 0.4 ^ Pf ^ 0.95. In this latter region, which is relevant for physics 
searches related to top jets for instance, the difference between the two approximations is 
significant. In this region, the strongly ordered approximation fails to capture much of the 
essential physics. This has implications for parton-shower calculations of this quantity: 
we do not expect unmatched parton-shower calculations to be accurate. Matching to tree- 
level matrix elements — so long as it is done to a sufficiently high multiplicity as to ensure 
at least three matched partons inside a jet — will introduce the required corrections to 
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the strongly ordered splitting functions. However, even the matched calculations may be 
quite sensitive to the matching procedure, and in particular to the size of the remaining 
region where the pure shower calculation is used. 

D. Behavior at Large Planar Flow 

Coming back to the 1 — > 3 computation (53), let us consider the jet function at Pf = 1. 
At this point 

A = - m^elel sinV (m^ + 2^2^102 cos (pf < , (61) 

and we are within the integration domain only when A = 0: the phase space dimension 
is reduced by 1. Because the splitting function has no singularities within this domain,^ 
it implies that the jet function vanishes at Pf = 1 in our approximation. 

Fig. 3 confirms this by showing that the semi-analytic jet function drops to zero as the 
planar flow approaches 1. As we will see below, the drop is a feature of our three-body 
approximation, and it will not be present when higher order corrections in as are included. 
It also shows the fit to the leading- log result, including the sub leading coefficients in 
Eq. (48), and we notice that the semi-analytic jet function diverges from the fit at the 
level of 10% near Pf = 0.95. We heuristically take this point to mark the beginning of the 
drop. Our three-body approximation is not valid beyond this point. 



^ All splitting-function singularities come from soft or coUinear limits, which imply Pf — >• 0. 
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• Quark jet 
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FIG. 3: Semi-analytic jet function near Pf = 1. The dashed hne shows the fit to the leading 
form of Eq. (48) with the coefficients given in Sect. VB. 
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E. Comparison of Running Scales 
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FIG. 4: Ratios of gluon jet functions with different choices for fi. The jet functions are evaluated 
at p = 1 TeV, m = 180 GeV and R = 0.4, and they are divided by the jet function that uses the 
hybrid (dipole) scale. 



In this section we consider how different choices of the running scale /i affect the jet 
function. Recall that the jet function (30) is of C(a^); ^^^^ evaluate the two powers 
of as at different scales. One power we evaluate at the jet-mass scale (corresponding to 
the first 1 — )■ 2 splitting in case of hierarchical emissions); but there are several possible 
choices of scale for the second power, corresponding to the second (softer) 1 — )• 2 splitting. 
We consider three possibilities: 

1. Set /i to be the jet mass m, 

2. Set /i = minjjjsy}, and 

3. Set yU to be a hybrid scale, described by the 2 — )■ 3 dipole scale for gluon emission 
or Sij in the g ^ qq case; see appendix E for details. 
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We expect the last of these to be the most accurate one, and we compare the others against 
it in Fig. 4. The choice of scale makes a significant difference to the value of planar-flow 
jet function, of order 10-30% at relevant values of the planar flow. Furthermore, the 
ratios are not constant as a function of the planar flow. A variation is to be expected in a 
leading-order calculation, as nothing in the matrix element compensates for the change of 
scale. We would expect this variation to be substantially smaller in a next-to-leading order 
calculation of the jet function. The scale variation is hidden in parton-shower calculations 
(as each algorithm chooses one particular scale), but this should be considered as an 
intrinsic source of uncertainty. Unlike the error made by applying a strongly ordered 
approximation, this uncertainty is not removed by matching to tree-level matrix elements. 
Matching to one- loop matrix elements as well would be required to reduce it. 

The results displayed in Fig. 4 at large Pf can be understood in a simple way. At Pf = 1 
we have a symmetric configuration of partons, where Sij = iv? /?>. Let us assume this is 
the dominant configuration. For a three-gluon configuration (which dominates the gluon 
jet function), the dipole scale is then /Xdipoie = ^/3. We expect that the jet function ratios 
at Pf = 1 will be given by the corresponding ratios of couplings. 



as(min(sjj)) as{m 



2^ 



(/^dipole) ' ' (/^dipole; 

This agrees nicely with Fig. 4, within numerical uncertainties. 



0.92, "\ ' . ^ 0.84. (62) 



VI. COMPARISON TO PARTON-SHOWER CODES 



In this section we compare our semi-analytic result to parton-shower simulations of QCD 
scattering at the LHC. The simulations are of 2 — )■ 4 (matched) and 2 — 2 (unmatched) 
matrix-element scattering, plus showering. We used MadGraph/MadEvent 5 [ ] with 
Pythia 6.420 (virtual ordered) showering [ ], and SHERPA 1.3.1 [ 1 1], and in both cases 
we have used the CTEQ6L set for the parton distribution function [ ) ]. The jet algorithm 
is anti-/cT [16] with R = 0.4, implemented in Fast Jet [17, 48]. The other parameters 
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are ^/s = 7 TeV, 950 < px < 1050 GeV and \ri\ < 1. We integrate over the mass 
window 160 < m < 200 GeV, which is consistent with our motivation for this work. The 
simulations include showering but not hadronization or detector simulation. As we show 
below in Sect. VIA, the effect of hadronization at large planar flow is to increase the 
distribution by about 15%, while leaving the shape of the distribution unchanged. 

The semi-analytic results are computed at p = 1 TeV and integrated numerically over 
the same mass window. There are separate jet functions for gluon and quark jets, and 
the total jet function is given by J = xJg + (1 — x)Jq, where x can be thought of as the 
"fraction of gluon jets" in the sample. While this is not a well-defined quantity, we may 
get a rough estimate for it by considering matrix-element 2 — )■ 2 scattering. Then x is 
given by the ratio of outgoing gluons to total outgoing partons (within the and r] cuts), 
and using this method we find x ~ 0.24. 

The parton-shower distributions are normalized such that the integral of each over the full 
jet-mass and planar-fiow ranges is 1. In particular, this means that the area under the 
plots presented below is not 1. Our jet functions are naturally normalized in the same 
way, so we expect the total jet function to agree with the parton-shower result within its 
range of validity. 

Fig. 5 shows the comparison of the semi-analytic results to the parton-shower jet functions 
for MadGraph (with Pythia showering) and for SHERPA. The second factor of the strong 
coupling as is evaluated at the hybrid scale as described in Sect. VE. Away from the 
peak we find that the parton-shower results fall between the semi-analytic quark and 
gluon functions, in agreement with the theoretical prediction. Notice that the parton- 
shower jet functions display no special behavior near Pf = 1, while our semi-analytic jet 
functions drop to zero there (see Sect. V D). This discrepancy is due to missing higher-order 
contributions in the theoretical calculation. In detailed comparisons with our theoretical 
result we will exclude this highest-Pf region, restricting ourselves to the range Pf < 0.95. 
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FIG. 5: Jet functions from parton-shower simulations, and our predictions of subsection VB 
using 1 — )• 3 splitting functions for the quark and gluon planar-flow jet functions. The upper 
(green) solid curve shows the gluon jet function, and the lower (black) solid curve the quark one. 
The dashed (green) curve is the average jet function with gluon fraction x = 0.24. The region of 
expected validity of the semi-analytic form is highlighted. The points with error bars show the 
parton-shower results, in the highlighted region from top to bottom: matched SHERPA (red), 
unmatched SHERPA (yellow) , matched MadGraph (dark blue) , and unmatched MadGraph (light 
blue). 

The simulations include an infrared cutoff of 15 GeV at the matrix-element level, which 
represents the minimal k-r distance between pairs of partons.^ A SHERPA simulation 
with a higher cutoff of 25 GeV (not shown in Fig. 5) gives a qualitatively similar result. 

One can check, for example by generating random three-body jets, that a cutoff of 15 GeV 
implies the matrix-element results apply only at Pf > 0.4. Below this only 1 — )■ 2 splittings 



^ In SHERPA, the CKKW matching scale Qcut also serves as this cutoff. 
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are in effect, and in addition we expect that resummation and non-perturbative effects 
become important. In support of this, below we will show that the effect of hadronization 
in the simulation does not alter the shape above Pf ~ 0.4. On the theory side, as we 
approach the peak (at Pf ~ 0.1) from above, resummation effects become important and 
our perturbative approximation breaks down. For the purpose of comparison we will 
therefore restrict ourselves to Pf > 0.4 . 

Fig. 6 shows a detailed comparison of the parton-shower and theoretical jet functions. 
The region in which we expect to find agreement is highlighted. Finally, Fig. 7 shows the 
comparison to the theoretical result using the iterated 1 — )■ 2 splittings (see Sect. VC), 
with the second factor of the strong coupling as evaluated at the scale of the second 
splitting Sij. Comparing with Fig. 5, it is clear that using the 1 — i- 3 splitting function 
results in a significantly better approximation to the jet function. (The first factors of 
the strong coupling are evaluated at the jet-mass scale in both cases, but the choice of 
second scale is different.) 
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FIG. 6: Detailed comparison between jet functions from parton-shower simulations and our 
semi-analytic results. The figure shows (J — Jq)/{Jg — Jq) where Jg and Jq are respectively the 
semi-analytic gluon and quark jet functions. These jet functions are shown by the upper (green) 
solid line at 1 and the lower (black) solid line at 0. The average jet function with gluon fraction 
X = 0.24 is shown by the dashed (green) line. The region of expected validity of the semi-analytic 
form is highlighted. The points with error bars show the parton-shower results, with the order 
and color coding as in Fig. 5 . 
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FIG. 7: Comparison of jet functions from parton-shower simulations and semi-analytic results 
using a strong-ordering approximation (iterated 1 — )• 2 splitting functions). The two factors of 
the strong coupling Og are evaluated at the jet-mass scale and the scale of the second splitting 
Sij. The region of expected validity of the semi-analytic form is highlighted. The curves and 
points are as in Fig 5 . 



A. Hadronization 

The effect of turning on hadronization in the parton-shower simulation is shown in Fig. 8. 
Below Pf = 0.4, where matrix-element events are discarded due to the infrared cutoff, we 
see that hadronization affects the shape of the jet function significantly. Above this value, 
hadronization affects only the overall normalization. 
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FIG. 8: A comparison of parton-shower jet functions with and without hadronization, shown as 
a ratio of the two. Both distributions are from MadGraph + Pythia simulations with matched 
jets. 



B. Estimation of Non-Collinear Corrections 

Our theoretical computation relies on two approximations. The first is working to leading 
order in perturbation theory, which is valid for sizeable mass and planar flow, well above 
the peak locations in the resummed distributions. The validity of this approximation 
requires that resummation effects be negligible, as{m) \og{m/pR) ^ 1 and as{m) log Pf ^ 
1, for m and Pf respectively. It also requires that the values of these variables not be too 
close to kinematic boundaries, such as m = pR and Pf = 1. The second approximation 
is the collinear approximation, valid for narrow jets. The two approximations limit the 
region of validity, while offering a non-trivial window of applicability for the calculations. 
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The window is home to a variety of potential new-physics searches, for example those 
using boosted top-quark jets. 

We would like to have a better understanding of the corrections to our theoretical results, 
and accordingly we would like to separate the collinear corrections from those due to re- 
summation. For that purpose we consider a matrix-element calculation, with no showering 
or hadronization, and compare it with our theoretical prediction (see Fig. 9). This calcu- 
lation is carried out at leading order in a^, and thus differs from our theoretical prediction 
only in corrections to the collinear approximation, as the tree-level matrix elements em- 
ployed are exact throughout phase space. The difference between the two results provides 
an estimation of the corrections away from the i? — )■ limit. We see in Fig. 9 that these 
corrections vary significantly with Pf. 

We normalize the exact leading-order jet function to match the way the semi-analytic 
jet function is normalized at leading order ((9(a^)), so that 

-1 

dpT^drfdm? dPi 

Here, o"2-s>/,x is the 2 — )■ 2 cross section (x denotes a parton of any flavor), and a2^j,x is 
the 2 — )• 4 cross section for an outgoing state that includes a jet J of matching type. Both 
cross sections are computed to leading order in a^. The integrals are over the window 
950 < Pt < 1050 GeV, |r/| < 1 and 160 < m < 200 GeV, with the jet function itself 
evaluated at = 1 TeV. 

Let us relate these objects to quantities that are directly measurable in a Monte-Carlo 
integration. Consider the differential cross section da in a given bin, for example the cross 
section for events within our kinematic window, and with planar flow in a small range 
[Pf,Pf +rfPf). It is given by 

dN , , 

da = —a, (64) 

where dN is the number of events inside the bin, is the total number of events produced 
in the integration, and a is the total cross section computed by the simulation. Using this 



/dm^J^f^ = ( [ dp^dri ^'^'^^f'^ [ f dpTdridm^ - — '^^^T^'^'^ ■ (63) 
^ \J dp^dr] J J -^--^-'-^^t 
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relation in Eq. (63), we find that the jet function is given by 

/ dm J/(Pf) = — ^ . (65) 

Here, N2-^x and (y2^x are the total number of events and total cross section produced 
in the 2 — )■ X simulation; A^2-5-/,a; is the total number of events in the 2 — )■ 2 integration, 
with a parton of flavor /, that fall within our kinematic window; and N2^j,x is the total 
number of events in the 2 — )■ 4 integration, with a jet of flavor /, that fall within our 
kinematic window (including jet mass), and within our Pf bin. 




+ 



+ 



4- + 



MadGiaph^ matrix -element quark jets 



Pf 

(a) (b) 

FIG. 9: Comparison of gluon and quark semi-analytic jet functions, with the results from Mad- 
Graph matrix-element 2 — )• 4 calculations. The parameters are p = 1 TeV, m = 180 GeV and 
R = OA. 



VII. CONCLUSIONS 

The planar-flow (Pf) distribution of highly boosted narrow massive jets is interesting, 
because a non-vanishing value of this variable implies that the corresponding jet consists 
of at least three hard partons in a perturbative description. QCD jets with sizeable planar 
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flow and jet mass (m) form an important background to various new-physics signals. 
For instance, massive jets with large Pf arise in models with heavy resonances decaying 
dominantly to top quarks or in supersymmetric models with i?-parity- violating gluinos. In 
this paper, we have studied the planar flow distribution of narrow QCD jets. We obtained 
a semi-analytic form for this distribution, independent of the underlying hard process 
giving rise to the jet. We have made use of QCD factorization properties to do so, and 
have computed jet functions which express the probability of a parent parton fragmenting 
into a jet of given planar flow and mass. We computed the leading-order approximation 
to these jet functions using the universal 1 — > 3 tree-level coUinear splitting functions. We 
compared this approximation to a strongly ordered coUinear approximation, using iterated 
1 — > 2 splitting functions, and flnd substantial differences. Our results are, unsurprisingly, 
sensitive to the choice of scales in the strong coupling. We have also derived the leading- 
log behavior of the jet functions analytically. Our results arc expected to be valid only in 
the range 0.4 ^ Pf ^ 0.95, and for sizeable jet mass, as fixed-order predictions will diverge 
both as Pf ^ and m — >■ 0. The divergence of the planar-flow distribution should be 
regulated by resummation of leading logarithms to all orders in perturbation theory. To 
the best of our knowledge, this resummation has not been computed, and the resulting 
resummed distribution would be of interest. 

We have compared our semi-analytic jet function to parton-shower predictions using var- 
ious event generators. The broad features are in agreement in the region of validity of 
our flxed-order calculation, as the parton-shower results interpolate between the predicted 
quark and gluon jet functions. The details differ, however, between the parton-shower 
results and a suitably-weighted average of quark and gluon jet functions. We find that the 
results from SHERPA are above our semi-analytic calculations, both with and without 
matching (using CKKW). In contrast, the results from MadGraph with Pythia showering 
are above our results with matching (using MLM), but below our predictions without 
matching. In the region of validity, hadronization effects increase the value of the jet 
function modestly without altering its shape. We note that the results of these two 
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parton-shower codes do not agree with each other, either with or without matching. The 
differences between them, even with matching, are comparable to the differences from our 
semi-analytic results, and the differences are greater without matching. This variation 
suggests that some caution should be exercised when comparing these results with ex- 
perimental data, and that a data-driven approach to jet substructure should be explored 
as well. The qualitative agreement of the semi-analytic results with the parton-shower 
calculations does however suggest that futher refinement of the fixed-order prediction, 
for example by carrying out a next-to-leading order calculation, would be valuable. The 
required 2 — )■ 4 one-loop matrix elements are available, and have already been used for 
phenomenological studies [49, 50]. 

Planar flow is one of several three-prong substructure variables that can play a role in 
discriminating highly-boosted jets arising from three-body decays of heavy particles from 
QCD backgrounds. It tends to be sensitive to soft radiation near the edge of the jet, 
which also makes it sensitive to pile-up. This motivates the use of jet filtering or template 
overlaps, where the hard substructure in a jet can be enhanced in a controlled manner. 
Such enhancement would be expected to bring jet functions closer to the fixed-order per- 
turbative ones calculated in this paper. 
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Appendix A: 1 — 2 Splitting Functions 



In this section, for completeness, we quote the 1 — )■ 2 spin-averaged and color-averaged 
splitting functions [33, 38] P/1/2 at leading order in a^. Let us define 



The functions are 
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Appendix B: 1 — 3 Splitting Functions 



In this section we write down the 1 — )■ 3 spin-averaged and color-averaged splitting func- 
tions -P/1/2/3 at leading order in a^. We follow the conventions of Ref. [iS]. Our definition 
of Zi agrees with this reference in the coUinear limit. In addition to Eq. (Al), we define 
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Appendix C: Three-Parton Kinematics in the Narrow- Jet Approximation 



In this appendix we write down kinematic quantities of three-parton configurations, at 
leading order in the coUinear approximation 9i < R <^ 1, m terms of the integration 
variables of Eq. (33). As dictated by the kinematics we assume m < pR (see Eq. (35)), 
and this implies that m <^ p. 

For parton three-momenta we use spherical coordinates {pi, 6i, (pi) with i = 1,2, 3, relative 
to the jet axis p, and we define = 0i — 02- In the coUinear approximation, 
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where 6ij is the angle between partons i and j. 
For parton 3, we have 
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In expanding the square root in we assumed that p — pi — P2 is large relative to the 
second term. This assumption in not valid in the limit of soft ^3, where in fact p — Pi— P2 
becomes arbitrarily small (and even negative). To avoid this limit we restrict ourselves to 
the range pi < P2,P3- 

The angles 6ij are given by 
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Finally, for the jet mass and planar flow we have 
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Appendix D: Analytic Leading-Log Coefficients 



In this section we write down the coefficients ^/ of the leading-log jet function (48) 
computed analytically from the integral (47). We find 
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The constants C^, Cp and Tr are defined in Eq. (Al). 



49 



Appendix E: Strong Coupling Renormalization Scales 



In this section we explain in detail our choice of running scale jj for the second factor 
of as that appears in the jet function (33). The first factor, corresponding to the first 
1 — 7- 2 splitting, is evaluated at the jet-mass scale. To make a realistic choice for fi (which 
corresponds to the scale of the second 1—7-2 splitting) we use the dipole model [52-56], 
in which a 2 — )■ 3 splitting of partons is described as an emission from a color dipole 
consisting of the two parent partons. The natural scale for this process is given by^° 



in the case where partons 1 and 3 form the dipole, and parton 2 is being emitted. In our 
case we do not know which of the partons is emitted in the second splitting, and in fact 
the second splitting is not even well-defined in general. We choose the scale to be the 
minimal one among the several options (corresponding to permutations of the partons), 
relying on the splitting function's preference for soft-coUinear emissions. For example, in 
the case oi g ^ ggg splitting, the scale is given by 



The dipole model does not describe cases where the second splitting process is quark-pair 
production, and in such cases we choose fi to be the mass of the produced pair. The scales 
for the remaining processes are given by 
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In the limit of soft-collinear emission, fi becomes the pT of the emitted parton relative to its parent. 
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Appendix F: Mass and Planar-Flow Leading-Order Jet Functions from Splitting 

Functions 

In this appendix we consider the jet function's behavior in the hmit of small planar flow, 
basically using an iterated 1 — > 2 splitting analysis in the limit of strong ordering. Let 
us first examine in more detail how one can obtain the jet-mass distribution in the soft- 
coUinear limit: the single emission rate is proportional to, 

, ^, dz dO , 
doi^2{z,9) ^ —— (Fl) 

where z is the energy fraction of the emitted gluon and 9 is the angle between the emitted 
gluon and the parent parton. In this approximation the jet mass is 

m] = E'^z9'^ (F2) 

For a fixed emission angle, 9, we can rewrite Eq. (Fl) as 

da{m],9) oc F3 
im?j 9 

where now > 9^ > m?j/E'^. We see that for a fixed mass, the distribution of 9 is 
characterized by 1/9. The jet-mass distribution is obtained upon integrating 9 between 
the boundaries of integration, giving 

dai^2 1 , fR^E'^ 



*„5 

where in hadronic collisions E should be replaced by px- The missing proportionality 
coefficient is nothing but as Ca/t^ (for a gluon jet). 

Next let us try to obtain an expression for the planar-flow distribution in the limit of small 
planar flow for a massive jet. In this limit, we expect that the dominant contribution arises 
from configurations where the third parton is soft and collinear to either of the first two 
partons. This contribution should be described well by a 1 — )■ 2 splitting function. We 
can thus iterate the above expression starting with a two parton configuration of mass mj 
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where the two patrons arc separated by an angle ~ d. We can now add a second emission 
with an angle 9' such that 9' <^9 and an energy fraction z' <^ z such that the third parton 
can be thought as soft and coUinear with respect to either of the first two partons. The 
differential cross section describing this configuration is given by, 

1 / 2 n\ 1 /ni i\ dmjd9 dz' d9' 

da{m% 9) x dai^2{9', z') ex -^-^ ^ (^5) 

where now z' is integrated between and z/2 at most. (At the upper boundary, both 
partons have the same energy; this violates the soft approximation but is highly suppressed 
due to the weight in the splitting function.) Likewise, 9' is integrated between and 9 at 
most. For a fixed value of mass and planar flow z' and 9' are not really independent to 
leading order. The tensor I describing this configuration is 

^ ^^f [z, (1 - z)]z'9"sm'<j> [z9s, (1 - z){l - 9.,)]z'9' sm2(p/2\ 

\[z9,, (1 - z){l - 9s)]z'9' sm2<p/2 (1 - z){9 - 9sf + z9l ) ' 

It is easy to check that = (1 - z)9 and (1 - z){9 - 9sf + z9l = z{l - z)9^ = mV-E^ 
and the bracketed expressions in the (11,12,21) entries cover the cases where the third 
parton is emitted from the harder or softer of the first two patrons, the softer one being 
characterized by an angle = (1 — z)9. (f) describes the azimuthal angle — the third 
parton's emission angle relative to the line connecting the two hard partons. To leading 
order in z' we find 

Pf = 4— [z, (1 - z)]z9'^ sin^ ^ 4— z'9''^ sin^ , (F7) 

where on the right-hand side of the above relation we have focused on the most singular 
region z' <^ z <^ 1. Similarly to the jet-mass case we can interchange z' and Pf in the 
singular region, Pf <^ 1. Integrating over the first emission and using the splitting function 
for both radiations (assuming only gluons for simplicity), focusing on the most singular 
region z' <^ z <^ \, 9' <^ 9 ^ \ and changing variables from z' to Pf (we can ignore as 
in this approximation the splitting function does not depend on it) we find 
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For a given Pf the range for 9' is 

I — 2Tn 

\/Pf— (F9) 

Integrating over 9' yields the final expression for the Pf distribution for jets of small Pf 
and mass, 

/ (fa \ alCl^ fR^E^\ I f ER \ 

fejp,,, - ^^^^ l^j ^ Pf i2WPf J 

where in hadronic collisions E should be replaced with px- 
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